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Abstract 



This paper surveys a few aspects of the global theory of wave equa- 
tions. This material is structured around the contents of a minicourse 
given by the second author during the CMI/ETH Summer School on 
evolution equations during the Summer of 2008. 

1 Introduction 

The week-long minicourse on which this brief survey paper is based came 
after a vigorous, detailed and outstanding series of lectures by Jared Wun- 
sch on the applications of microlocal analysis to the study of linear wave 
equations. Both lecture series took place at the Clay Mathematics Institute 
Summer School at ETH Zurich in 2008. The goal of this minicourse was to 
describe a few topics which involve global aspects of wave theory, relying at 
least to some extent on the microlocal underpinnings from Wunsch's lectures. 
The first of these topics is an account of some striking consequences that 
can be derived from the finite propagation speed property. While this had 
been applied in various interesting ways before, the systematic development 
of this principle appears in the very influential paper of Cheeger, Gromov 
and Taylor |CGT82j . We recall how this property, applied to solutions of the 
wave equation associated to a Laplace-type operator, can be used to obtain 
estimates for solutions of various related operators. We present only one 
application of this, which is a lovely argument due to Gilles Carron which 
estimates the off-diagonal decay profile of the Green function for generalized 
Laplace-type operators on globally symmetric spaces of noncompact type. 
This result had caught the lecturer's eye in the months before this Clay 
meeting and nicely illustrates the unexpected power of the finite propaga- 
tion speed method. Following this, the remainder of the lectures reviewed 
several diifferent approaches to scattering theory and described a few of the 
relationships between these. The primary goal, however, was to introduce 



the Friedlander radiation fields and explain how they give a concrete real- 
ization of the Lax-Phillips translation representation. We follow suit here, 
recalling the outlines of a few of the numerous successful approaches to 
scattering theory and culminating in a discussion of these radiation fields. 

This paper attempts to give some feel for what was presented in these 
lectures. The reader should be warned that the topics covered here are in 
many places old-fashioned and we omit any mention of many of the most 
important recent advances and trends in scattering theory. The material 
here is meant to indicate a few things that can be accomplished, often with 
not very sophisticated machinery by modern standards. We typically make 
very restrictive assumptions in order to convey the main essence of the 
ideas. We give references for further reading interspersed inter alia, but do 
not make any claim to a comprehensive bibliography. 

The material assembled here is based on the notes of the first-named 
author; the lecturer (and second author) is extremely grateful to him not 
only for this careful recording of the lectures, but also for his enthusiasm 
during the lectures and his very substantial assistance in writing this paper. 
We did discuss at some point, but later abandoned, the possibility of writing 
a much more exhaustive treatment of some of the topics here, particularly 
the theory of radiation fields. That will unfortunately have to wait for 
another day and other authors. We hope that this survey accomplishes 
what the original lectures also attempted, which is to whet the reader's 
curiosity to learn more about this subject. Needless to say, wave theory is 
an immense subject and we mention here only a very small set of possible 
topics. 

Throughout this paper we focus on properties of solutions, and of the 
solution operator, for the wave operator 

Uv = D'^-L, where L = V*V + 1/ (1.1) 

acting on sections of some bundle E over a Riemannian manifold {M,g), 
where V is the covariant derivative of some connection on E and y is a 
(self-adjoint) potential of order 0, which can either be scalar or an endo- 
morphism of E. For simplicity we typically assume that V is smooth and 
compactly supported, although neither of these properties are present in al- 
most any of the interesting physical or geometric applications. Furthermore, 
we often discuss only the scalar Laplacian and its perturbations, although 
the extension of all results below to this slightly more general framework is 
usually just notational. Finally, here and below we write D = jd. 

As noted above, we take advantage of the luxury of being able to refer 
back to the excellent lecture notes by Jared Wunsch |Wun08j covering his 
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longer minicourse. Those notes provide a nice introduction for many central 
themes and results in the subject, including the existence of solutions of 
the equation Dn = / with vanishing Cauchy data, or of Dn = with 
prescribed nonzero Cauchy data, along a noncharacteristic hypersurface, 
the positive commutator method leading to Hormander's renowned theorem 
on propagation of singularities of solutions, the finite propagation speed 
property, and much else besides. Using this as a blanket resource, we can 
dive right into the material at hand. 

There are now many terrific monographs concerning the local and global 
aspects of wave equations. Michael Taylor's three-volume series |Tayll| 
belongs high on this list; it contains an amazing amount of information 
about many different topics. Other recent monographs with a particular 
focus on hyperbolic equations include those by Alinhac |Ali09| , Lax |Lax06] , 
Ranch [Ran 12] : we mention also the new book by Zworski on semiclassical 
analysis |Zwol2j . 

The authors are very grateful to the Clay Foundation for making this 
Summer School possible - it was a lot of fun and the large attendance 
and enthusiasm of the participants was amazing. We also appreciate the 
forebearance by the editors of this volume for their (relative) tolerance for 
the length of time between the original lectures and when this paper was 
finally written. Both authors are very grateful to many people for teaching 
us about many of the topics here. We thank, in particular, Gilles Carron, 
Richard Melrose, Gunther Uhlmann, Andras Vasy and Jared Wunsch. Gilles 
Carron and Andras Vasy also gave some helpful remarks on this paper. 

D.B. is supported by NSF Postdoctoral Fellowship DMS-1103436; R.M. 
is supported by the NSF grant DMS-1 105050. 

2 Finite propagation speed and its consequences 

Although jWun08| contains a proof of the basic finite propagation speed 
property for the operator Dy, we begin by recalling this familiar argument 
very briefly. We then show how using the functional calculus one can write 
the Schwartz kernels of various functions of the elliptic operator L in terms 
of the Schwartz kernel of the wave operator. This leads directly to the im- 
portant Cheeger-Gromov- Taylor theory which uses finite propagation speed 
to obtain interesting estimates for these Schwartz kernels. We illustrate this 
with an outline of Carron's estimates for the resolvent and heat kernel of 
generalized Laplacians on symmetric spaces of noncompact type. 
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2.1 Finite propagation speed 

The fundamental identity behind finite propagation speed is the observation 
that for any sufficiently regular function u, 

div^iutVu) = utOou + ^dtiu^t + |Vup). (2.1) 

We suppose that the space on which we are doing calculations has a global 
time function t and moreover, splits as MxM, with a static Lorentzian metric 
—dt^+h, where (M, h) is a Riemannian manifold. A hypersurface Y C MxM 
is called spacelike if its unit normal v (with respect to this Lorentzian metric) 
satisfies v ■ u < Suppose that $7 C M x M is a domain bounded by two 
spacelike hypersurfaces, dO. = Yi U which meet transversely along a 
codimension two submanifold, and that u is a solution of the homogeneous 
wave equation, Dqu = 0. Integrate (12. Ih over Q,. The left side is transformed 
using the divergence theorem; the first term on the right vanishes while the 
second term is also transformed to a boundary integral. If Vj = {vt,j, ^x,j) is 
the upward-pointing unit normal to Yj, decomposed into its vertical {t) and 
horizontal {x) components, then we obtain 

/ (|utp + |Vup)|i/t,i| - 2ut ■ du^u\vx^i\ 




Since Vj is timelike, the integrand on each side is bounded from below by 
c(|utp + |Vnp) for some c > which depends on Yj. We conclude that 
if ut = Vn = on Yi, then these same quantities must also vanish on 
Y2- Finally, if is foliated by spacelike hypersurfaces, then the vanishing 
of (utjVu) on the bottom (spacelike) boundary of 0, can be propagated 
throughout this entire region, and hence if u vanishes at Yi, then li = in 

n. 

If we consider wave operators with terms of order or 1, then this calcu- 
lation can be adapted to show that if O is foliated by spacelike hypersurfaces 
Zg, < s < 1, then the integral over Zs of |Mtp-|- |Vnp satisfies a differential 
inequality, and the fact that it vanishes when s = implies that it vanishes 
for all s < 1. 

To interpret this calculation, we observe that there many natural do- 
mains which can be foliated by spacelike hypersurfaces in this way. Indeed, 
suppose that p = (ti,xi) is any point, and P^" denotes the (backward) 
domain of dependence of this point, i.e. the set of points in M x M which 
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can be reached by timelike paths traveUing backward in t and emanating 
from Let Y be one of the level sets {t = to} where to < ti- Then 

the region Q, = {{t,x) £ T^Ti xi '■ ^ — ^0} can be shown to have a space- 
like foliation by submanifolds Yg which all intersect along the submanifold 
{{t, x) £ xi ■ ^ ~ ^o}- Thus any homogeneous solution of Dyu = which 
vanishes along with its normal derivative along {t = to} vanishes through- 
out this Q. This implies that if the Cauchy data of u at to is supported in 
some subset K, then the Cauchy data of n at ti = to + where r > 0, 
is supported in the subset K^- = {(ti,x) : distg{x, K) < r}, which is pre- 
cisely what is meant by saying that the support of a solution propagates 
with speed 1. For more general variable coefficient hyperbolic equations, 
the speed of propagation may be variable but is still finite. 

2.2 Cheeger-Gromov- Taylor theory 

Consider the fundamental solution for the problem 

Dyu = 0, u\^^Q = <j), dtu\^^Q = 0. 

It is customary to write this solution operator as cos(t\/L), so that the 
solution u{t,x) is equal to cos(t\/Z)(/). We assume for simplicity that L has 
no negative eigenvalues so there is no ambiguity in this expression. This is 
an instance of the functional calculus for self-adjoint operators, which are 
defined in purely abstract terms using the spectral theorem, on one hand, 
but which can be used to describe solution operators for various equations 
involving L. There are many interesting examples, including prominently 
the resolvent and heat operator 

Rl{X) := {L - X^y^ and e"*^. 

The abstract definitions of these operators (i.e. defined using the spec- 
tral theorem) are all well and good, but in order to use them one usually 
wishes to know much more about their mapping properties. For example, a 
priori, using only these abstract definitions, we only know how one of these 
functions of L acts on functions, but not on other function spaces. The 
goal then is to obtain a more concrete understanding of the Schwartz kernels 
of any one of these operators. Of course, there is a lot of theory devoted 
to doing just this. Thus the classical theory of pseudodifferential operators 
gives a nice picture of the resolvent for any one value of A, while the theory 
of semiclassical pseudodifferential operators provides a means to understand 
this family of operators as A tends to infinity in various directions in the com- 
plex plane. Similarly, the well-known heat-kernel parametrix construction. 
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cf. |BGV92j ■ gives a way to understand the asymptotic behaviour of the 
Schwartz kernel of the solution operator for the heat equation in various 
regimes of the space M"^ x M x M. These theories and constructions give 
very precise information, but are often very intricate, and furthermore, it 
is often hard to use these ideas directly to say anything interesting about 
global behaviour of these Schwartz kernels. The idea in |CGT82j is that 
one can extract, often in a rather simple way, some very useful global be- 
haviour of these kernels using mainly the finite propagation speed property 
of cos(i\/L) and some other simple properties, such as the fact that the 
norm of cos{ty/L) as a bounded operator on never exceeds 1. 

To explain this, suppose that f{s) is a smooth, even function on M which 
decays sufficiently rapidly so that the following manipulations are justified. 
Assuming L > for simplicity, we define f{VL) using the spectral theorem, 
but at the same time we can spectrally synthesize this function of L directly 
from the wave kernel: 



1 /•oo 

/(VI) = — / f{s)cosisVL)ds. 



The simple but crucial observation is that this is not just an identity about 
abstract self-adjoint operators, but also calculates the Schwartz kernel of 
/(\/L) in terms of the Schwartz kernel of the wave operator. 

The following discussion is drawn from the paper [CGT82j . Suppose 
that / has the property that its Fourier transform /(s) is integrable, along 
with a certain number of its derivatives, on M \ (— e, e) for any e > 0. The 
first key result is that under such a hypothesis, if it E has support in a 
ball Br{y), then for R> r, 

\\f{VL)u\\L2(^M\Bu{y)) <'^'^\\u\\l2 / \f{s)\ds. 

J R-r 

The proof is very simple. We know that cos(s\/Z)u has support in B^j^\g\{y), 
so that 

||/(\/L)n|| < - / f{s)cos{sVL)uds < -||n|| / \f{s)\ds. 

A very similar argument gives bounds for \\LP f{y/L)L'^u\\ depending on 
the integral of some higher derivatives of /(s) over the same half-line. The 
particularly useful aspect of this is that the integrals of which appear on 
the right in these estimates start at R—r rather than at 0, and hence if these 
functions decay at some rate, then the right sides of these inequalities exhibit 
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the corresponding decay. Assuming we are on a space with appropriate local 
uniformity of the metric (or coefficients of L) , then we can deduce from this 
some off-diagonal pointwise estimates for the Schwartz kernel f{^/L){z,w). 
By off-diagonal we mean that the estimates are valid in any region where 
dist {z,w) ^ 0. One reason for assuming this local uniformity for L is that 
these arguments require bounds on the injectivity radius and volumes of 
geodesic balls, for example, in order to pass from to pointwise estimates. 

2.3 Carron's theorem 

This subsection provides a concrete example of how this all works. We 
describe some of the main features in the paper [CarlOj of Carron which uses 
the ideas above to derive fairly accurate pointwise bounds on the off-diagonal 
decay of the resolvent kernel and heat-kernel for Laplace-type operators on 
symmetric spaces of noncompact type. 

In order to describe this we must first explain at least a small amount 
about the geometry of these spaces. This is recounted elsewhere in much 
greater detail; the classic reference is |Hel84| . but we refer (self-servingly) to 
|MV05j for an analyst's point of view of this geometry. 

Symmetric spaces are distinguished amongst general Riemannian mani- 
folds by the richness of their isometry groups. Their defining property is that 
the geodesic reflection around any point {expp{v) i— )■ expp(— u)) extends to 
a global isometry; Cartan's classic characterization is that any such space is 
necessarily of the form G/K, where G is a semisimple Lie group and K C G 
is a maximal compact subgroup, endowed with an invariant metric. Because 
of this, almost all of the basic structure theory can be reduced to algebra 
and hence described quite explicitly. We shall focus on one particular real- 
ization stemming from the polar decomposition G = KAK, where K is as 
above and A is a maximal connected abelian subgroup. For a symmetric 
space X of noncompact type, this subgroup A is isomorphic (and isometric) 
to a copy of R'' for some k, where the positive integer k is called the rank 
of X. Using this polar decomposition, we identify G/K = KA. The map 
<I> : X A — 7- X, <I?(/c, a) = ka, is surjective, but far from injective. 

It is best to think of the simplest special case, the real hyperbolic space 
H"; here A = R and K = SO(n). The image of the origin G ^ via $ is 
a single point o G X, and this point is fixed by the entire (left) action of 
K. The space X is the union of geodesic lines through o which all intersect 
pairwise only at this point. The group K acts transitively on this space 
of geodesic lines through o with stabilizer SO(n— 1). Note that there are 
elements of K which take a geodesic to itself but reverses its orientation; 
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this means that we get a less redundant 'parametrization' by restricting <I> 
to K X R+. Geometrically, we have the familiar picture of H" as M"'' x 
with the warped product metric dr"^ + sinh^ r dO"^ . 

For a general symmetric space X of rank /c > 1, this picture generalizes 
as follows. The space X is the union of the various images of A by elements 
k £ K, and all of these images intersect at o, though kAOk'A often consists 
of a larger subspace. These translates of A hy k £ K should be thought of 
as the radial directions in X. Another important piece of structure is the 
existence of a finite set of linear functionals A = {oj} on A called the roots. 
These divide into positive and negative roots, A = A"*"UA~, and the positive 
roots determine a (closed) sector V C A hy V = {aj > y aj £ A+}. This 
sector V is the analogue of the half- line in H", and the restriction of to 
K X V is still surjective, and ii K' C K denotes the isotropy group at a 
generic point, then we can regard X as being the product V x [K/K') with 
certain submanifolds of K/K' collapsed along various boundary faces of V . 
In terms of this data, we can finally write down the multiply-warped product 
metric 

g = da^ + sinh^ aj dUj , 

j 

where the sum is over positive roots, da'^ is the Euclidean metric on A 
and drij is a metric on a certain subbundle of the tangent bundle of K/K' 
corresponding to the root aj. 

For simplicity here we just discuss the scalar Laplacian A on X and, 
following the theme of this section, consider the problem of estimating the 
Schwartz kernels of /(\/^^A) for suitable functions /. Because A commutes 
with all isometrics on X, the Schwartz kernel Kf{x,x') of this operator 
depends effectively on a smaller number of variables. Given any pair of 
distinct points x, x' € X, choose an isometry <^ of X so that (p{x') = o and 
(p(x) lies in some particular copy of A. If we ask that (p{x) = a £ V C 
A, then </? is almost uniquely determined. We thus have that Kf{x,x') = 
Kf{ip{x),o) = Kf{a,0). In other words, Kf is really only a function of the 
k Euclidean variables a = (oi, . . . , a^). In particular, when the rank of X is 
1, then Kf reduces to a function of one variable r > 0. 

This reduction points to the difficulty of studying functions of the Lapla- 
cian on symmetric spaces of rank greater than 1. Indeed, while the resolvent 
kernel R{X; x, x') on a space of rank 1 depends only on dist(x, x') and hence 
can be analyzed completely by ODE methods, the same is not true when 
the rank of X is larger. Similarly, even in the rank 1 setting, the heat ker- 
nel H{t,x,x') depends on two variables, t and dist(x,x'), but unlike the 
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Euclidean case, there is no extra homogeneity which reduces this further 
to a function of one variable. Thus the problem which Carron's theorem 
answers is how to give good estimates on these reduced functions, R{\, a) 
and H{t, a), where a € A is the 'relative position' between x, x' S X. 

Theorem 2.2 (Carron [CarlO] ). Let X be a symmetric space of noncompact 
type and rank k and consider the Schwartz kernels R{X, a) and H(t, a) of 
the resolvent (—A — Xq — A^)~^ and heat operator e*^, written in reduced 
form as above. The number Aq here is the bottom of the spectrum of —A; it 
may be calculated explicitly. Then 



The function '5t(a) is a somewhat messy but quite explicit and understand- 
able function which is a rational function of a and certain powers of t. The 
linear functional p on A is half the sum of the 'restricted ' positive roots; this 
is a standard object which appears frequently. 

It is known that the upper bounds given here are sharp in the sense 
that there are lower bounds that differ just by the constant multiple for 
these same kernels. We refer also to the papers |AJ99j and |LM10j sharper 
bounds obtained by different and more complicated methods. 

The proofs of these estimates are clever but not very long, and in the 
remainder of this section we give a few of the ideas which go into them. 

The first step is that if a G A is arbitrary and e G (0, 1), then we can 
estimate from above and below the volume of the set KB (a, e), where B{a, e) 
is a ball of radius e in A centered around a. The reason this can be done 
is because we have very good information on the Jacobian determinant for 
the coordinate change implicit in some natural coordinatizations induced by 
K xA-^ X. 



Let us first study the resolvent. Fix a € A such that dist(a,o) > 2. 
Let D be the "annular" region KB{a, 1) and suppose that a G L'^{D) (so a 
vanishes outside D). We estimate R{X,a)a and then R{X,a) itself, outside 
this region, i.e. where dist(a, o) > 2. If u = R{X, ■)a, then 



and 
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Then using that || cos(^-v/— ^ — Xo)\\l'2^l'2 = 1 as well as finite propagation 
speed, because of the support properties of a, we obtain that 

poo -Rc(A)5 1 



Jdist(a,o)-2 \^\ 1^1 

From here, using local elliptic estimates, we obtain that 

\u{o)\ < C7e-^°W'^'"*("'°)||(j||i2p). 

In other words, this estimates the norm of the mapping T defined by [D] B 
a R{X)a\^, whence (using the — ^ L°° norm of TT*), 

\R{X,x,o)\'^dx < C7e-2R°W^'^*('^'''). (2.3) 



D 



We next wish to find a similar estimate where the integral on the left is 
only over some ball B(ka, 1/4:) C D rather than the entire annular region 
D. More specifically we assert that 

Vol(S(A;a,l/4)) [ \R{X,x,o)\'^ dx < Ce-^^<^'>'^'''^'''°\ 

JB{ka,l/4) 

This must be true since if it failed for all such balls then the sum over all 
balls would lead to a violation of (j2.3p . 

Finally, noting that the volume of this ball is approximately e^^^"'^ , and 
applying the same local elliptic estimates as before to estimate the value at 
a point in terms of a local norm, we conclude that 

\R{\,a,o)\ < Ce~i^^W^*'^*("'°)-^'('^). 

This is the desired off-diagonal decay estimate. 

The corresponding argument to estimate the off-diagonal behaviour of 
the heat kernel proceeds in a very similar way, substituting local parabolic 
estimates for local elliptic estimates. We refer to jCarlOj for details. 

It is worth remarking that there are other very effective ways to estab- 
lish so-called Gaussian bounds for heat kernels under rather general cir- 
cumstances. We mention in particular the beautiful theory developed by 
Grigor'yan and Saloff-Coste, see |SC02j . |Gri09j . These techniques work in 
far more general circumstances, and depend on quite different underlying 
principles. However, one point of interest in Carron's work is that he is able 
to obtain the correct polynomial factors in front of the exponential, which 
is perhaps impossible using those more general approaches. 
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3 Scattering theory 



For the second and longer part of this survey, we turn to an entirely different 
aspect of the global theory of wave equations and discuss some approaches 
to mathematical scattering theory. This classical subject has deep physical 
origins, and has received numerous mathematical formulations. While these 
approaches are mostly equivalent, the correspondences between them are 
sometimes not always obvious. In the following pages we first review one 
point of view on stationary scattering theory, then turn to some perspec- 
tives on the corresponding time-dependent theory. This is all done with a 
distinctly PDE (rather than, say, operator-theoretic) focus. We conclude 
with a discussion of a more abstract functional analytic setup of scattering 
theory due to Lax and Phillips centered around the notion of a translation 
representation and explain how the theory of radiation fields developed by 
Friedlander provides a concrete realization of the translation representation. 

There are numerous settings in which to introduce any of these top- 
ics, including scattering by potentials, which is the study of Schrodinger 
operators —A + V on M", or scattering by obstacles, which studies these 
same operators but on exterior domains M" \ O with some elliptic boundary 
condition at dO. There are also significant differences between these theo- 
ries when the ambient space is odd- or even-dimensional. Finally, it is also 
natural to consider these same problems on manifolds which are asymptot- 
ically Euclidean or asymptotically conic at infinity (or indeed, have some 
other type of asymptotically regular geometry, e.g. asymptotically hyper- 
bolic). Each setting requires different sets of techniques, and in order to 
make this exposition as simple as possible, we focus on the combination of 
hypotheses where everything works out most simply. Namely, we study the 
scattering theory associated to L = — A + y on an orfd-dimensional Eu- 
clidean space M", with the strong assumption that V £ C^. We describe 
the structure theory for solutions of the Helmholtz equation (L — A^)n = 0, 
and for Dy :=□ + !/ = — L, the time-dependent wave equation, and 
give some indication how objects in these respective settings correspond to 
one another. 

There are very many excellent references to each part of what we discuss 
(and much that is closely related that we do not discuss), so we relegate 
almost all of the technicalities to those sources. We mention in particular 
[RS781 Vol. IV], ITayllj Ch. 9], |Per83| . |YaflOj and |Mel95] . The material 
on radiation fields is spread over several papers, starting from the original 
work by E.G. Friedlander jFriSOj . There is a forthcoming and detailed survey 
of this subject by Melrose and Wang |MW| . to which the discussion here is 
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intended to be an introduction. 

3.1 Stationary scattering theory 

The stationary formulation of scattering theory concerns the elhptic operator 
L — A^, where here and below, L = —A + V, with V G Cq° (and real- valued!). 
It is obvious that L is bounded below, i.e. 

/ {Lu)udV >-C I \u\^dV 

for all u € C^(M"'), and with little more work one can also prove that it 
has a unique self-adjoint extension as an unbounded operator on L^(]R"). 
Indeed, this is yet another consequence of the finite speed of propagation, 
see |Che73j . Its spectrum is contained in a half-line [— C, oo); the positive 
ray [0, oo) comprises the entire continuous spectrum, and there are a finite 
number of eigenvalues in the [— C, 0). If we allow V to be less regular, 
simple examples show that this negative interval may contain an infinite 
sequence of such eigenvalues converging to 0; the basic example of this is 
when V{x) = —l/\x\, which is the potential for the Schrodinger operator 
modeling the hydrogen atom. 

Assume initially that A lies in the lower half-plane QA < 0. Provided 
that A 7^ — «|Aj| corresponding to any of the negative eigenvalues —Xj < 0, 
the operator L — has an bounded inverse, 

Ry{X) = {L-X^)-\ 

This is called the resolvent and is a meromorphic family of bounded opera- 
tors on with simple poles at the points — i|Aj|. 

The first issue is to show that the continuous spectrum (A^ G [0,oo)) is 
absolutely continuous, or in other words, that the singular continuous part 
of the spectrum is empty. More specifically, we must find an L-invariant or- 
thogonal splitting L2(M") = Tipp Ti ac, SO that the restriction of L to ^^pp 
is discrete, while the restriction of L to H^c 

is absolutely continuous. It is a 
classical theorem due to Friedrichs that in this setting any eigenvalue of 
L is strictly negative. The proof consists of showing that if any such eigen- 
value is positive, then the corresponding eigenfunction must vanish outside 
a compact set, which violates standard unique continuation theorems. (This 
uses that V is compactly supported - ifV only decays rapidly then the argu- 
ment is a bit more intricate.) By the general spectral theorem, the absolute 
continuity of is equivalent to the existence of a unitary isomorphism 
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U : Hac — > L'^{^]Y), where Y is an auxiliary Hilbert space, so that the 
self-adjoint operator U o L o on L^(]R;y) is multiplication by the co- 
ordinate function t G M. One of the goals of scattering theory is to exhibit 
this unitary isomorphism explicitly, which is done using the M0ller wave 
operators, see below. A closely related goal is to understand the structure 
of generalized (non-L^) solutions to the equation {L — }?)u = 0, > 0. The 
key tool for all these questions is the resolvent i?y(A), introduced above. 

Let us first consider the free Laplacian Lq = —A on M". When A G 
1R\ {0}, the nullspace £{X) of the operator —A — A^ contains the plane wave 
solutions e*'^^ '^ for any oj G S"^^. Any linear combination of these plane 
waves also lies in ^(A), and indeed, general superpositions of these plane 
wave solutions span all of £{\). We explain this more carefully. For any 
5 gC°°(S"-1), define 

u{z)= [ e'^'-'^g{co)duj. 

This is a solution of (—A — A^)it = 0, and the most general (polynomially 
bounded) element of £{X) can always be obtained from this same represen- 
tation but allowing 51 to be a distribution. The "smooth" elements of <f (A) 
are those where g is smooth. 

We can look at this a different way. Note that since w 1— )• z • w is a Morse 
function on S"~^, and has critical points u = ±z/\z\, the stationary phase 
lemma shows that (assuming g is smooth), the integral expression for u has 
an asymptotic expansion of the form 

00 00 
u{z) ~ e^^l-'llzl""^ l^r^a+jW + e-'^^'^\z\-^ (3-1) 

j=0 j=0 

Here z = \z\6, 6 G S"~^ are polar coordinates on M". As part of this, one 
obtains that up to a multiple of 27r, a±^o = i^^^~^^/'^g{±6). Closely related 
is the assertion that any u G <f (A) has an expansion of this same form and 
moreover, fixing any a+^o ^ P'(5'"~^), there is a unique u G £{X) with this 
distribution as its leading coefficient. It is reasonable to regard the operator 
V : a+^o ^ u as solving a Dirichlet problem at infinity for —A — A^, and 
hence we call V the Poisson operator. 

The free scattering operator at energy A is the map 5o(A) sending the 
function a+^o to a_^o- Using the explicit representation above, we see that 
in this free setting, So{X)a{9) = i'^~^a{—9); it is just a constant multiple of 
the antipodal map. 
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Proceeding slightly further with the free problem, suppose that Im A < 0. 
Using the Fourier transform, one can determine the inverse of —A — (as 
an operator on Schwartz functions) via 

Mm = (-A - A')"V = (2vr)-" / e-<(|Cp - X^y'fiO dC 

When n is odd, this has a particularly simple form: there is a simple poly- 
nomial Pn{oi) of degree (n— l)/2 such that the Schwartz kernel of i?o(A) can 
be written as 

k-^?~>n(A|z-z'|)e-^^l^-^'l. (3.2) 

(In particular, P3(a) is simply a constant.) There is a related but slightly 
more complicated formula when n is even. This explicit expression shows 
that as a function of A, -Ro(A) continues holomorphically from the lower 
"physical" half-plane {9A < 0} to the entire complex plane when n > 3 
is odd. When n = 1, this continuation has a simple pole at A = 0, and 
when n is even, there is a similar continuation but to the infinitely sheeted 
logarithmic Riemann surface branched at the origin. To make sense of this, 
one can say that this Schwartz kernel continues as a holomorphic function 
taking values in distributions; an alternate and equivalent sense is to regard 
the continuation taking values in the space of bounded operators — )■ L^^^, 
(this domain space consists of compactly supported functions). From 
(j3.ip and stationary phase, one proves that if / G C^(M"'), then 

i2o(A)/ = e-*^l^l|z|-Vu;, 

where w is a smooth function on the radial compactification of M". This 
last assertion about smoothness on the compactification is simply a concise 
way of stating that w has an asymptotic expansion 

oo 

w ~ Wj{6)\z\~^ . 

j=0 

Let us now pass to the analogous considerations for the operator L. Some 
versions of all the structural results about solutions remain true. These are 
typically proved by a perturbative argument, which means that one no longer 
has explicit formula. The starting point is the Lipman-Schwinger formula, 
which gives a relationship between -Ro(A) and i?y(A) in the region in the 
A-plane where they both make sense. This states that 

Rv{X) = Ro{X) {I + VRo{X)r^ = (/ + Ro{X)vr' Ro{X). 
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The issue is to prove that the inverses of / + VRo{X) and / + Ro{X)V make 
sense, and to do this one observes that VRq{X) and Ro{X)V are compact 
operators (between suitable function spaces), so that one can invoke the 
analytic Fredholm theorem to obtain that these inverses, and hence Rv{X) 
itself, are meromorphic on the region where -Ro(^) is holomorphic (hence on 
C when n is odd and greater than 1). 

The argument sketched earlier that L has no eigenvalues embedded 
in the continuous spectrum implies that RviX) has no poles on the real axis. 
(The argument for regularity at A = requires slightly more care.) On the 
other hand, the negative eigenvalues Xj of L correspond to poles of Rv{X) 
at — i|Aj|. The new and perhaps unexpected phenomenon in that i?y(A) 
may have poles in the upper half-plane (and indeed, this always occurs if 
V is nontrivial). These poles are known as the resonances of L, and their 
location and statistical distribution has been the target of much research. 

Let £v{X) denote the nullspace of L — A^ (say in S'{M."')). Just as in the 
free case, this space may be generated using "distorted" plane waves. These 
are defined as follows. For any uj G S"~^ and A S M\ {0}, there is a function 
W\^i^ which is smooth on the radial compactification of M" so that 

lies in SviX). Note that the second term here is simply RQ{X){—Ve^^^''^). 
Superpositions of these can be used as before to generate all elements of 
SvW- Indeed, if 5 S C°°{S^~^), then the general "smooth" element of 
<?y(A) can be written as 




Using stationary phase as before, this integral has an asymptotic expansion 
of exactly the same form as (|3.ip . The leading coefficient a_|_^o(^) is again 
just (a multiple of) g, but now the other leading coefficient a-fl{9) is not 
simply the reflection g{—9), but rather a sum of this reflection plus an extra 
term which is an integral over involving both g and V. The scattering 
operator 5y(A), which sends a+^o ^ f^-.Oi is again unitary, and is the sum 
of the antipodal operator and another term which has a smooth Schwartz 
kernel. The map VviX) which sends a+^o to u is again called the Poisson 
operator. 

The results and definitions above continue to hold in suitably modified 
form not only for obstacle scattering, but also in the rather general setting 
of asymptotically Euclidean or asymptotically conic manifolds (these are 
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called scattering manifolds [Mel94j by Melrose). For more on this as well as 
many further details about everything discussed above, we refer to the book 
of Melrose |Mel95| . see also |Mel94j and |MZ96| . 

3.2 Time-dependent scattering 

We now turn our attention to the time-dependent formulation of scattering 
theory, and its relationship with stationary scattering. This time-dependent 
theory involves the study of "large time" properties of solutions of the wave 
equation. The connection with the stationary approach is via the Fourier 
transform in time; indeed, this Fourier transform carries L-L»| to L- A^, 
and asymptotic properties of as |i| — )• oo correspond to 'local in A' properties 
of the latter operator. For the wave equation associated to L = —A -|- V, 
where V is compactly supported, the intuitive picture is that one sends in 
a wave for times t <C from some direction at infinity and then observes 
what happens as this wave interacts with the potential and then scatters 
into a sum of plane waves as t +00. Amongst the many good sources 
for this material, we refer to the books of Friedlander [Fri75j . Lax [Lax06j . 
Lax-Phillips |LP89j . Taylor |Tayll| and Melrose |Mel95j . 

3.2.1 Progressing wave solutions 

We begin by describing the special class of progressing wave solutions for 
wave operators. The calculations here go back to the dawn of microlocal 
analysis and can be regarded as the nexus of many constructions and ideas 
in that field. This construction is quite geometric and it is most naturally 
phrased in terms of the wave operator on a general Lorentzian metric g. 
The special case of a static metric g = — dt^ + h on the product of M with 
a Riemannian manifold (M, h) is of particular interest, and we discuss at 
the end how this specializes for the particular operator Dy = □ -|- y on 
Minkowski space. 

Thus let (X, g) be a Lorentzian manifold and consider -|- V ^ where 
V G C^{X). We look for solutions u to (Dg -|- V)u = which have the form 

u = ip a(r), 

where 93 is smooth, a is a distribution on M which models the 'wave form' 
of the solution, and F is a function on X with nowhere vanishing gradient 
which we call the phase function. To be concrete, we typically let a = 5 or 
a = for some A; > 0, but the key feature we require of a is that it behave 
like a homogeneous function in the sense that its successive derivatives and 
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integrals are progressively more or less smooth than a itself. Of course, 
it is usually impossible to choose solutions of (D^ + V)u = which have 
this precise form, but the goal is to add increasingly higher order correction 
terms of a similar form involving the integrals of a so that, in the end, this 
initial expression is the first term in some asymptotic expansion of an exact 
solution. 

The first step is to calculate 

(□3 + V)u = -^di (g'^V\i\dju) + Vu 
Vm ^ ' 

= a"{r)g (Vr, Vr) ^ + a'{r) {2g (Vr, Vip) + ipDgF) 

+ a(T){Dip + Vip) 

As indicated above, assume that is a sequence of distrubutions on M 
such that Ok = a^+i- (Again, refer to the basic example oq = 6, Ok+i = 
^x^.) Let us now assume that 

u r^'^Uk = '^(pkit,z)ak(T). (3.3) 

fe>0 fe>0 

We apply the calculation above and group together the terms of the same 
order (where the order of is k and each derivative lowers the order by 1). 

Grouping terms of the same order, we attempt to choose ip^ so that each 
term vanishes. The only term of order —2 is (poaQ{T)g (VF, VF), so the first 
requirement is that 

5(vr,vr) = 0. 

This is known as the eikonal equation and states that VF is a null-vector 
for the metric g. This is a global nonlinear Hamilton-Jacobi equation for F. 
In the special case X = M x M, g = — dt^ + h, the eikonal equation can be 
written as 

if we write F = t — S", where S is a function on M, then 

|V,,S|2 = 1. 

It is straightforward to see that the level sets S = const are at constant 
distance from one another, so in general, S{x) = dist/^ (x, Z) where Z is 
some fixed level set of S. Even in the more general Lorentzian setting, the 
function F incorporates a lot of the distance geometry of g. 
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In any case, fix a solution T of the eikonal equation. We liave now ar- 
ranged that the term of order —2 vanishes. In fact, for any k, the term 
containing 5'(Vr, VF) vanishes, and so the equations for the higher coeffi- 
cients simplify to transport equations. In particular, the term of order — 1 
reduces to 

Since VF is nowhere vanishing, this is a linear ODE for ipo along the integral 
curves of VF, which means that given any initial conditions for (pQ on the 
characteristic surface F = constant we may solve this equation locally. 

The term of order k — 1 yields an inhomogeneous transport equation for 
ipk in terms ofT,ipQ, . . . ,{pk. We solve this transport equation with vanishing 
initial data and proceed inductively to choose all fk- 

It is possible to asymptotically sum the series (j3.3p . This means that we 
can choose a function v with the property that 

N 

V - y^V9fcafc(F) 
fc=o 

is as smooth as the next term in the series, (Pn+ioln+i{^)- By construction, 
(□g + V)v = / € C°°{X). We must now invoke a theorem guaranteeing the 
existence of a smooth solution w for the initial value problem {^g + V)'w = f 
with vanishing Cauchy data vanishes, where / is smooth. Given this, then 
u = V — w is a solution of the original equation and the expansion we 
have calculated determines the singularity profile of u. Note that these 
singularities of u occur precisely along the union of level sets F = c where 
one (and hence every) is singular at c. 

For the special case where g = —dt^ + dx^ on Minkowski space, fix 
oj € S"~^ and consider the equation 

{d'j - y) u = 0, u = 5{t - z ■ oj) when t < 0. 

The eikonal equation |VF|^ = has solution F(t, z) = t — z • uj. This gives a 
global solution of the wave equation for all t when V = 0. However, by the 
propagation of singularities theorem, the wave front set of the solution u for 
the perturbed problem with this initial data in the distant past agrees with 
that of this exact free solution. Hence it makes sense to look for a solution 
of the perturbed problem of the form 

u ^ 6{t — z ■ uj) + ^k{t, z)x^{t — z ■ u), 

k>0 
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for some choice of smooth functions ipk- This fits exactly into the scheme 
above (and was, of course, the setting for the original version of these cal- 
culations). The first transport equation is 

2{dt-uj- V,) = 0, 

which means that 999 is a function oi t = z ■ oj and z] its Cauchy data is 
defined on the hypersurface t = z-uj^ and the equation dictates that it must 
be constant along the lines parallel to llj. 

Once we have determined lpq, . . . ,(pii., then the (A; + 1)^* transport equa- 
tion is 

2{k + i){dt-uj- V,) ifk+i = -{a + v) ifk, 

which we solve with vanishing initial data. Carrying this procedure out for 
all k determines the Taylor series of u along the hypersurface {t = z-co}. As 
described earlier, we can use the Borel Lemma to choose an asymptotic sum 
V for this series, so that {O -\- V) v = f is smooth and v satisfies the correct 
"initial condition" for t <^ 0. We can then find a smooth correction term w 
which solves away this error term. Thus u = v — w is a,n exact solution 

The calculations here were historical precursors to the more elaborate 
but ultimately very similar ones which come up in the construction of Fourier 
integral operators. Indeed, solving the eikonal equation for T is the direct 
analogue of solving the eikonal equation for the phase of an FIO. For po- 
tential scattering, keeping track of the parametric dependence on uj fixes 
the phase; the solutions of the transport equations are the coefficients in 
the expansion of the amplitude, and these correspond to the terms in the 
expansion for the symbol of the FIO. 

3.2.2 M0ller wave operators 

We now turn to another perspective on time-dependent scattering, which 
is through the definition of the so-called M0ller wave operators. This can 
be regarded as a formalization of the discussion above; there we described 
how to calculate the profile of the solution obtained by "sending in" a delta 
function along a particular direction. Our goal now is to put this information 
together into a map which compares the long-time evolution with respect to 
the perturbed equation against that for the free equation. 

Let us suppose now that g = —dt^ + h \s a. static Lorentzian metric. For 
any (C^) potential V , define the wave evolution operator 

Uv{t) : C~(M") X C^{W) C7~(M") x C^{W), 
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where, if u solves the Cauchy problem 

{a + v)u = o, {u,dtu)\t=o = {(l),i^), 

then C/y(to)(</', V') = iu,dtu)\t=to- The free wave evolution operator Uo{t) 
is defined analogously using solutions for Du = instead. Uniqueness of 
solutions of these Cauchy problems implies that Uy and Uq are groups, i.e. 
C/*(t)"^ = U,{-t) and U^t + s) = [/*(t)[/*(s) for * = or V. 
Now define the M0ller wave operators W± by 

W±{<t),^)= \\m C/y(-t)C/o(t)(</.,V), 

t— >±oo 

presuming the limit exists. This limit is meant to be taken in the sense of 
strong operator convergence. If we define the energy space 

He = ^{^,^I^): y"^' + |V,,/.|2dy <oo|, 

then W± extends by continuity to all of He- It can be proved that if certain 
local measurements of this energy decay appropriately, then —/S. + V has 
no eigenvalues and this extension is an isomorphism of He to itself. If 
— A + y does have eigenvalues, then T-Lpp determines a finite dimensional 
subspace in He and W± is an isomorphism from He onto the orthogonal 
complement of T^pp, which we denote H^. 

Since Uo{t) and Uv{t) are unitary, the wave operators W± are charac- 
terized by the property that 

\\Uv{t)W±{4>,i^) - Uo{t){(t>,i:)\\H^ ^ as t ^ ±00. 

for all V') £ H-^. Now define the scattering operator 

S = W^^W-\ 

this is an isomorphism of H-^. It describes the relationship between the 
asymptotic free wave emerging as t +00 for a solution of the perturbed 
equation (□ + V)u = in terms of the incoming free wave for t <S 0. 

These operators lead directly to the unitary isomorphism mentioned ear- 
lier which intertwines L (or rather, its restriction to Tiac), with a simple 
multiplication operator. In other words, the existence and properties of the 
wave operators and scattering matrix proves that the singular continuous 
spectrum of L is empty. 

There are many other settings where one can define analogues of the 
M0ller wave and scattering operators. Classically this is done for exterior 
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domains, and more recently on asymptotically Euclidean or conic manifolds 
(where the structure of the scattering matrix is quite intriguing, see |MZ96j ) , 
as well as other geometric settings such as asymptotically hyperbolic man- 
ifolds, etc. There is also a parallel but vigorous line of research concerning 
the possibility of defining the analogues of wave and scattering operators for 
various classes of nonlinear evolution equations. 

3.2.3 Lax Phillips theory and radiation fields 

hi this final section we present yet another approach to scattering theory. 
This is the more abstract approach developed by Lax and Phillips |LP89j . 
which has played an influential paradigmatic role. Directly following this we 
describe the theory pioneered by Friedlander [FriSO] on what he called the 
radiation fields associated to solutions of a linear wave equation. These de- 
scribe certain asymptotic information about waves, and beyond their purely 
analytic appeal, they also provide a beautiful realization of the Lax-Phillips 
theory. These radiation fields have received quite a lot of attention in re- 
cent years, and the theory has been extended to various nonlinear settings 
as well. There is a forthcoming and much more detailed survey specifically 
about radiation fields |MWj to which we direct the reader. 

Throughout this section we fix a Hilbert space % and a unitary semigroup 
U{t) which acts on it. The specific application we have in mind is that % is 
the space He of finite energy initial data for the wave equation on with 
n odd and U{t) is the wave evolution operator. More precisely, let T-Lq be 
the completion of the space C^(M") x C^(M"') with respect to the norm 

\UM?H, = [ (\v<p{z)f + mz)\') dz- 

then, for ((/>, V') G let Uo{t){<j),ip) be as defined in the previous section. 
The unitarity of Uq corresponds to conservation of energy for solutions of 
this wave equation. 

Return now to the general formulation. 

Definition 3.4. A closed subspace T> d % is called outgoing, respectively 
incoming, if 

(i) Uit)!) CV for t > 0, respectively t < 0, 

(ii) f]tmUit)'D = {0}, and 

(in) \j-jj{tyD = n. 
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In the example above, the space P+ consists of the pairs {(p, ip) G T-Lq 
for which the solution u{t,z) vanishes for \z\ < t when t > 0. Continuous 
dependence of solutions of the wave equation on initial data shows that P+ 
is a closed subspace. The first and second properties follow from the obser- 
vation that if {(pjip) S Tio, then by finite propagation speed, the solution of 
the wave equation with initial data U{s){(l),'ip) vanishes for \z\ < t + s. 

The fourth property is more subtle. For the unperturbed wave equation 
in odd dimensions, it is a consequence of Huygens' principle; in even dimen- 
sions, one may prove it using local energy decay, but it can also be proved 
fairly explicitly via the Radon transform. We say more about this later. 

The fundamental result of Lax-Phillips theory is the existence of a trans- 
lation representation: 

Theorem 3.5 (Lax-Phillips jLP89j ). Let U{t) be a group of unitary oper- 
ators on %, and T) an outgoing subspace with respect to U{t). Then there 
exists a Hilbert space fC and an isometric isomorphism 



such that = L2 ((0,oo);/C) and ^ o U{t) =Tto^, where {T{t)f){s) = 

f{s — t) is the standard translation action 0/ M on L^(]R;/C). The isomor- 
phism $ is unique up to an isomorphism of fC. 

The isomorphism given here is called an outgoing translation representa- 
tion of U{t). There is an essentially identical result giving an isomorphism 
<!>' which maps an incoming subspace P_ to ((—00, 0); /C) and intertwines 
U{t) with T[t). This is called an incoming translation representation. The 
auxilliary Hilbert space fC may be taken to be the same as for the outgoing 
translation representation, but of course the map is diff'erent than 

Returning again to the unperturbed wave equation in R", n odd, there 
is an explicit way to obtain the translation representations using the Radon 
transform. 

Definition 3.6. For any f G C^(M"), define the Radon transform 



where da{z) is surface measure on the hyperplane {z,9) = s. Clearly Rf G 



$ : H ^ L2((-oo,oo);/C) 
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A key property of the Radon transform for our purposes is that it is 
invertible and in fact the inversion formula is quite exphcit: 

where \Ds\ is defined by conjugating multiphcation by \a\ with respect to 
the Fourier transform. A remarkable fact, which can be proved by direct 
computation, is that R intertwines the Laplacians on M" and M, 

RAf = dlRf. 

We now define the Lax-Phillips transform: for n odd, and ((/>, V') S 
C^{W) X C^{W), let 

LP(0, = (^^y(^ {Rct>) {s,9) - (i?^) {s,e)) . 

Theorem 3.7. For n odd, the Lax-Phillips transform LP extends to a uni- 
tary isomorphism 

LF-.Ho^L^ (R;L2(S""1)) , 
and is a translation representation, 

(LP Uo{t)i(P, ^|^)) {s, 9) = {Tt LP(</), V)) (s, 6) = (LP(</., V)) {s - t, 6). 

One consequence of Theorem 13. 71 is that T-Iq splits as an orthogonal direct 
sum of the incoming and outgoing subspaces: 

no = V+®V-. (3.8) 

In particular, in this special case, the outgoing and incoming isomorphisms 
$ and are equal. 

Now consider the wave equation with potential. As before, assume that n 
is odd and V G C^{W) is real-valued. Choose R so that supp V C B{0, R). 
Let U {t) be the group associated to the Cauchy problem 

au + Vu = 0, {u,dtu)\t=o = {^,^1^), (3.9) 

i.e. U(t){4),ip) = {u{t) , dtu{t)) . Since V does not depend on t, there is a 
conserved energy, 



\\{u{t),dtum\l= / (\dtu{t,z)\^ + \Vu{t,z)\^ + V{z)\u{t,z)\^) dz. 

(3.10) 
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The Hilbert space % is the set of pairs {(f), ip) for which this energy is finite. 
It is not hard to see, using the Sobolev inequahty, that H. and T-Lq consist 
of the same pairs of elements, although the norm is different. The energy 
extends to the bihnear pairing on %: 



this is anti-symmetric with respect to the pairing (|3.1ip . The wave group 
U (t) can be regarded instead as the solution operator for the system 



We now make a simphfying assumption that L = —A + V has no 
eigenvalues, or equivalently, that A has no such eigenvalues. Without this 
assumption, the results below require a projection off the finite dimensional 
space Tipp. We refer to |LP66j for more details about how to proceed with- 
out this assumption. The advantage of this assumption is that now the 
energy ()3.10p is positive definite. 

For this perturbed problem, we define the incoming and outgoing sub- 
spaces 'D±^R d % to consist of those elements {4>,il^) so that Uo{t){(j),tp) 
vanishes in |z| <t + R for t > 0, respectively \z\ < —t + R for t < 0. Thus, 
in terms of the free incoming and outgoing subspaces, P±,_r = Uo{-^R)V±. 
The verification that these satisfy all the correct properties relies on the 
following 

Lemma 3.12. // f = ((^, ^/^) G V+^r, then Uo{t)¥ = U{t)i for t > 0; the 
analogous statement holds for I^-,_r, when t < 0. 

We now use this to show that P+,_r is an outgoing subspace for U (t) onH. 
Indeed, by this lemma, the first two properties follow from the corresponding 
properties of 2?+. For the third property, suppose we know that for any 
compact subset K C and any solution u of p.9p . we have 




(3.11) 



Consider now the operator 





lim \\u{t)\\% j^ := 
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This is called local energy decay, and is known to be true in many circum- 
stances. Now consider the initial data i = {(pjip) £ Ti with f -L |J U{t)T>^^fi 
with respect to the pairing (j3.1ip . Thus U{t){ _L 'D+^r for any t, and in 
particular, U{t){ _L P+^j? with respect to the standard pairing on x L^. 
This shows that C/o(— i?)C/(t)f _L P+ with respect to the standard pairing, 
and hence Uo{-R)U{t)f G P_ and Uo{-2R)U{t)f G V^^r. 

Consider now v{s,z) = U{s)UQ{—2R)U(t)i. By Lemma [3.121 v{s,z) 
agrees with Uo{s)Uo{—2R)U{t){ for s < and thus vanishes for \z\ < —s + R 
for s < 0. 

Now we bring in the local energy decay. This implies that for any e > 0, if 
t is sufficiently large then ||C/(t)f||^ b{5R) ^ ^- -^'^'^ such t, finite propagation 
speed implies both 

\\Uo{-2R)Um\E,Bm < ^' r(-2^)t^(0f|U,B(3i?) < 

Because the two equations and the initial data agree outside B{R), using 
finite propagation speed again, we get that Uo{—2R)U{t){ = U{—2R)U{t){ 
for \z\ > 3R and hence 

\\Uo{-2R)U{t)¥ - U{-2R)U{t)i\\E < 2e. 

Because U{t) is unitary with respect to (j3.1ip . applying U(2R — t) to the 
difi^erence shows that 

\\U{2R - t)Uo{-2R)U{t)i - f||^ < 2e. 

Finally, since t is large, 2R - t < and so U{2R - t)Uo{-2R)U{t){ = 
Uq{2R - t)Uo{-2R)U{t)i by LemmaElSl This shows that in fact 

\\Uo{-t)U{t)f-i\\j,<2e. 

Because Uo{—R)U{t){ G P_, the first term here is an element of "D-^t-R and 
thus vanishes for \z\ <t — R. Taking t even larger gives 

l|f|li^ <2e, 

and therefore f = 0. This establishes the third property. 

Theorem 13.51 asserts the existence of incoming and outgoing translation 
representations for the incoming and outgoing subspaces T^-,r and T>^^r. 
We shall give a a concrete realization of these using the so-called radiation 
fields. 

Our next goal is to show that a particular quantitative rate of local 
energy decay implies that the local energy actually decays exponentially. 
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Theorem 3.13. Suppose that for each compact subset K C M" there is a 
function exit) which tends to as t ^ oo, such that if the Cauchy data u{0) 
have support in K, then 

\\um%K<cK{t)\\um\l- (3.14) 

Then there are positive constants C and a depending on K such that if u{0) 
is supported in K, then 

\Ht)\\E,K<Ce-^'MmE (3.15) 

for all t > 0. 

The proof uses the compactness properties of the Lax-Phihips semigroup 
Z{t), which we introduce now. If P±^fi. are the orthogonal projections onto 
the orthocomplements of T)±^ii, then Z{t) is given for t > by 

Z{t) = P+^RU{t)P_^R. 

The local energy decay hypothesis in the theorem statement implies that, 
for t large enough, Z[t) has norm bounded by 1/2, and repeated application 
of Z{t) leads to the exponential decay. 

We are now in a position to introduce the radiation field of a solution 
u to the perturbed wave equation. The idea is to identify initial data for u 
with a normalized limit of the solution along outgoing (or incoming) light 
rays. As before, we start with the definition of these radiation fields for the 
unperturbed operator. 

Suppose that u solves Don = with intial data (</>, ^). Introduce coor- 
dinates s = t — \z\ and x = l^;]"^; these parametrize the family of outgoing 
light rays and the position along them. Now define the auxiliary function 

: Ms X (0, oo)^ x S"~^ M, vis, x, 9) = x""^ u f s + i, -Q\ 

\ X X ) 

Here is simply z in polar coordinates. Since xP'gu is nondegenerate at 
X = 0, the function v extends smoothly across x = 0. We then define the 
forward radiation field operator Tty. 

Tz+{(t>,ii>){s,e) = dsv{s,Q,e) 

The derivative of v is included here to make TZj^ an isometric isomorphism: 

n+:no^ L^{R X s"-i) 
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Furthermore, the Minkowski metric is static, so TZ^ intertwines wave evolu- 
tion and translation in s: 

n+Uo{T){^, e) = 7^^.(</., - r, e) 

Now observe that if f G then T^+f vanishes when s > 0. This follows 
from the unitarity of the radiation field operator, and the inverse image of 
those functions in L^(R x S"~^) supported in the nonpositive half-cylinder 
form an outgoing subspace P+. 

P+ = {7^;V:/(s,^)=0for s>0} 

Indeed, this is a closed subspace; the first and second properties follow 
directly from the fact that is a translation representation, while the 
third property follows from the surjectivity of 7^+. One may also define P_ 
via the backward radiation field ; this encodes information from solutions 
in the limit as t \ — oo. For the free wave equation, TZ+ has an explicit 
expression in terms of the Radon transform, and this can be used to show 
that no = 'D+(SV_. 

For the perturbed equation the forward and backward radiation fields, 
C±, are defined in the same way. We can also define the scattering operator 
A using the radiation fields by 

A = C+CZ^ 

Thus A maps data at past null infinity into data at future null infinity. The 
relationship to the scattering operator S introduced in Section [3.2.21 is that 

s = n^^c+cz^TZ- = ni^Au^ 

The conjugation of A by the Fourier transform in s corresponds to the 
scattering matrix employed by Melrose in |Mel94j . 

The radiation field exists and is a unitary operator in a variety of geo- 
metric settings. On asymptotically Euclidean spaces, this is due to Fried- 
lander jFriSOl iFriOlj and Sa Barreto |SB03t ISBOSj : on asymptotically real 
and complex hyperbolic manifolds it was proved by Sa Barreto |SB05j . and 
Guillarmou and Sa Barreto |GSB08j . respectively. In the asymptotically Eu- 
clidean and real hyperbolic cases, Sa Barreto and Wunsch |SBW05j prove 
that it is a Fourier integral operator with canonical relation given by the 
sojourn relation, a close relative of the Busemann function in each of these 
geometric settings. The radiation field has also been defined in certain non- 
linear and non-static situations. In particular, the first author and Sa Bar- 
reto show |BSB12j that it exists and is norm-preserving for certain semilin- 
ear wave equations in M^, while Wang [WanlOt IWanl 1] defined the radiation 
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field for tlie Einstein equations on perturbations of Minkowski space when 
the spatial dimension is at least 4. Forthcoming work of the first author, 
Vasy, and Wunsch |BVWj analyzes the s — )• oo asymptotics of the radiation 
field on (typically non-static) perturbations of Minkowski space. 
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